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Abstract 

In this pedagogical note it is demonstrated how the numeric value of fine struc- 
ture constant may be established by comparing results following from the cal- 
culations in the framework of Quantum Electrodynamics with the experimental 
data. As an observable, the coefficient of 1/r in the Coulomb law is used. 
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Introduction 



It is a wide-spread perception that Quantum Field Theory [1], the framework for operat- 
ing with the fundamental degrees of freedom, serves as the basis of of our understanding of 
the Universe. According to Quantum field Theory a full information about the physical sys- 
tem under consideration is provided by the Lagrangian of the corresponding field-theoretical 
model. The only unknowns in this approach are the values of the parameters appearing in 
the Lagrangian, e.g. the electron mass and coupling constant appearing in the Lagrangian of 
Quantum Electrodynamics. 

In effective theories the values of parameters are assumed to be known. E.g. when analyzing 
the problem of the hydrogen atom spectrum using the Schrodinger equation, the electron and 
proton mass and the value of the fine structure constant are given; it is assumed that these 
values will be established in some "high-brow" , fundamental theory. Quantum Field Theory is 
the most promising candidate for this task. Of course, it is impossible to derive values of the 
parameters (like charge and mass) in the framework of Quantum Field Theory alone - some 
numeric values, serving as an input data arc necessary. Then, the best bet in establishing the 
values of the parameters would be the extraction of their values by comparing results following 
from the Quantum Field Theory with these input data. This procedure is described in the 
present note on the example of Quantum Electrodynamics the Lagrangian of which contains 
two parameters - fermion mass m and the coupling constant e. 

We will concentrate on e, assuming that the value of parameter m is known. Establishing 
the numeric value of m is a slightly easier task from theoretical viewpoint since mass is the 
invariant of Poincare group and does not require for its definition interacting fields, while 
couphng constant can not be introduced for a free fields - at least first non-trivial order of 
perturbation theory or some other approximation to an interacting quantum fields description 
is required. 

To proceed with the comparison of the results of Quantum Field Theory with the experi- 
mental data we need to make sure that the field-theoretical expressions contain no divergences. 
Mathematical method used to demonstrate that the results of calculations in Quantum Field 
Theory may be finite is called renormalization. Renormalization procedure is sketched below. 



As is well known, the available methods of calculation in the framework of Quantum Field 
Theory encounter infinities that either persist, leading to nonsensical results, or cancel among 
themselves, leaving finite parts [1]. In order to operate with the well defined quantities the 
regularization is introduced. The idea is to distort the original field-theoretical expressions so 
that mathematical manipulations become well defined. At the end of calculations the distortion 
parameter is set to the initial value. 

As an example, consider integral in the four dimensional space-time: 



In the regime \k\ 3> m we have J(l, 4) pa / d^k/k"^ which leads to an infinite result when 
integrating over all space. 



Brief description of renormalization 
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There are many ways to regularize divergences. In case of integral J(l, 4) one can introduce 
cut-off \k\ < A, or change the power of (A;^ — m^) from —1 to — s, or change the dimensionality 
of space time from 4 to any number n. Correspondingly, instead of J(l, 4) we will have: 



When A is finite, or when s > 2, or when n < 2 integrals do not diverge and the mathematical 
operations involving these integrals are well defined. After the desired result is expressed in 
terms of regularized quantities (in our case, in terms of J(4, 1; A), or J(4, s), or J(n, 1)), the 
"distortion parameter" is set back to its initial value: A = oo, or s = 1 or n = 4. 

In Quantum Electrodynamics (QED), containing only two input parameters - coupling 
constant e and the electron mass m, prescription leading to a finite results can be described as 
follows [2] : use regularization (say, dimensional regularization, prescribing to use n dimensional 
space-time at the intermediate stage of calculations [3]) to calculate any two physical quantities 
T,a and S^: 



where eo and mo are the input, unrenormalized coupling constant and the electron mass and 
the momentum dependence, being irrelevant for the discussion, is omitted. Standard choice for 
Eq, and Tip are the renormalized mass and coupling constant [1]: 



E„ = m = mo(l + e^zW + e^Z^^) +...), = e = eo(l + Z« + ef^Z^f +•■■), (2) 



but, in general, any two physical quantities can be used. 

In most cases, when the regularization is removed, expressions for and diverge. 
Therefore, regularization can not be removed, i.e. the limit n — > 4 does not exist in (1). Since 

and E/j are physical quantities, they must be finite^. The ansatz implicitly embedded in 
renormalization procedure is that and E^ are finite and it is just the imperfect math which 
prevents to demonstrate the finiteness of physical quantities. 

Theory is called renormalizable if there are finite number of physical quantities Ej, j = 
1, 2 • • • for which regularization can not be removed and the remaining Eat, Eat+i, • • • do not 
diverge after the regularization is removed. In other words, in renormalizable theory calculation 
of Ejv, Ejv+i, • • • leads to an explicitly finite expression. Theory is called non renormalizable if 
there is no physical quantity which remains finite after the regularization is removed [1]. 

Now, back to QED. It is renormalizable theory as it was demonstrated more than half a 
century ago [2]. 

Renormalization procedure for QED prescribes to solve equations (1) for eo and m.Q in terms 
of Eq, E/j and the regularization parameter n, and substitute into the expression for any other 



Definition of physical quantities of Quantum Field Theory is a highly non-trivial problem; e.g. 
contrary to what may seem self-evident, the matrix element of elastic electron-electron scattering 
diverges - only the matrix element accounting for the emission of infrared photons are finite, and thus 
can be interpreted as a physical quantities, representing a measurable [1]. 




= Ea(eo, mo; n), E^j = E;3(eo, mo; n) 



(1) 
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physical quantity S^. In QED expression for becomes finite when the regularization is 
removed, i.e. the hmit n ^ 4 exists: 

= hmE^(eo,mo;n) = limE-j,(eo(Sc,,E^;n),mo(Ilc,,E^;n);n) < oo (3) 

After the divergences are removed it becomes possible to compare the expression for a physical 
quantity with the numeric value for that quantity extracted from the experiment. This com- 
parison allows to establish the numeric value of renormalized parameters. E.g. the numeric 
value of renormalized coupling constant e is extracted according to following scheme: 

cT(m, e)th = (Texp e = fia^xp) (4) 

where Cexp is the measured value of physical quantity a and ath is the expression for a obtained 
in the framework of QED. This procedure is described below. 



Extracting value of e by comparing QED calculations with the Coulomb law 

As an observable a let us use the field created by a static source with the mass M ^ m 
where m is the electron mass. The proton at rest can be chosen as such a source. Prom the 
experiment it is known that at large distances, r ^ h/mc, the electromagnetic field generated 
by a proton is: 

Aoit, r) = ^, A(t, r) « (5) 

where 

„o 47r 

~ — (6 
^137 ^ ^ 

Equation (6) can not be derived from QED; using (6) as an input, we have to calculate the 
value of e from QED. 

We need to calculate the long distance asymptotic of electromagnetic field created by a 
static source, compare the result expressed in terms of e and m to the experimental information 
(5)- (6), and from that comparison extract the numeric value for e. 

Prom QED it follows that field created by the static source = eSfj,oS{r) is 

^(r) = J D'^''{r-r')j,{r')d^r' ^ e J dtD{v') ^ e J d^ke^^'' D{Q,k) (7) 

In (7) D = D^^ is a scalar part of the photon propagator D'^" and e is the QED coupling 
constant. 

In a first non-trivial order of perturbation theory expression for D is (from now on we omit 
inessential overall factors and use units h — 1, c — 1): 



/' 

Jo 



,2 



M ^ ~ IJ,^X(1 - X) 

dxx{l-x)[ln - + 
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+ 



/i^ (k^ + n'^)x{l -x) 
— ii'^x{l — x) 



(8) 



where 



a 



e 

An 



(9) 



and e = Co + 0(6^) is a renormalized coupling constant. 

(8) follows from the standard perturbation theory [1] and the straightforward algebraic 
manipulations; for our goal it is enough just to sketch the derivation and explain what is /i and 
how did it appear in (8). 

Calculation of the photon propagator requires the knowledge of the vacuum polarization 
operator 11^^,^: relation between these two is [1] 



9nv 



(10) 



Photon propagator is obtained by inverting equation (10). Polarization operator diverges; in 
order to obtain finite result for a photon propagator it is necessary and sufficient to get rid 
of infinities in n(A;^) where n(A;^) is a form-factor of the polarization operator: 'n.^y{k) = 
(^^ 9ij.u — k^k,^)Il{k'^). This may be achieved by subtracting twice at some normalization point 

k' 



1^ 



2. 



(11) 



That is how the (arbitrary) normalization scale ii appears in 11 and through (10) appears in 
the expression for the photon propagator (8). Exact expression (in case of perturbation theory 
- summed up series) for a physical quantity does not depend on an unphysical parameter // [1]; 
since wc use the first order approximation, (8) contains fi explicitly. 

Standard choice for the subtraction point in QED is /i^ = 0. Expression for the photon 
propagator appearing in a textbooks [1] is obtained using /x^ = 0; for our purposes we have 
to consider the case of arbitrary ^ and then find out which values of /i allow to calculate 
long-range asymptotic of in QED. 

Now, since we need asymptotic at large r, let us expand D in k^: 



' dxx{l-x) flnfl- ^'""^^ 



X 



+ 



/x^x(l — x) 



H^x{l 



+ 



X] 



+ 



11^ x{\ 



X] 



xil — x) 



k^ (777,2 _ j2'2x{l — x)) 

First two terms in 2q;[ J/tt of (12), 



/■I , ,^ , /, / it^x(l-x)\ iJ,'^x(l-x) \ 

/ dxx{l-x)[lnil- '^ \ M + 2 2 n J 
^0 \ V / — ijl^x[1 — x) J 



(12) 
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are k-independent; we denote the corresponding Fourier transform by 



2aF{m, /x) 
7rr 

Third term of 2a[ j/vr is ~ k^^; denote Fourier transform by 

li"^ ip(m, ii) r (14) 
Fourier transform of the fourth term, ~ k^, results into 

-^W (15) 

Coefficient in (15) is /i independent which means that the corresponding term should be directly 
interpreted as an observable. This is really the case: (15), known as the Uehling contact term 
[4] is observed in spectroscopical experiments. 
Collecting these together in (7) results in 



Ao{r) = e 



r 47r V 7rr loTrm^ tt / 



Q{a, m, iJ,, r) 



(16) 



Comparing the QED result (16) with the experimental information (5) demonstrates that 
Q{a, m, ^, r) is renormalization invariant and thus can be compared with Q. the coefficient of 
1/r in the Coulomb law. Note that so far the value of e is not known; we need to simplify (16) 
as much as possible and then compare the resulting expression to (5)-(6). 

Expression (16) contains term proportional to r - 2aii^(fi{m, /i) t/tt. Therefore, in calculating 
asymptotic for large r, terms 0{a^) can not be neglected: it becomes necessary to sum up all 



the powers of a, equate the result of the summation to y47r/137, and solve from this equation 
for the expansion parameter a = e'^/An. In other words, if the terms proportional to r are 
present in expression given by a fixed order of perturbation theory, it is necessary to sum 
up all the terms in all orders of perturbation theory and obtain expression for a physical 
parameter Q{a, m, /i, r) in a closed analytical form. Result will contain all the powers of a 
which means that the numeric value of the expansion parameter a can not be extracted from the 
numeric value of Q{a, m, /i, r) - it is a highly non-trivial task to solve a from Q{a, m, /i, r) — 
E°li a^Qjim, II, r). 

In QED there is a way out of this difficulty: the condition ^ 0^ allows to subtract at 
^2 — 0. When /x^ = the term ~ r drops out, and from QED it follows that in the limit 
r — > oo and when /i^ — field is 



^Examining expression for the photon propagator (8) demonstrates that D diverges when fj? = and 
when m? = 0. This problem of infrared divergences can be resolved in massless electrodynamics, but, 
it presents an imanswered challenge to Quantum Chromodynamics - theory of interacting quarks and 
massless gluons [6]. 
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^o(l")|,i2=o,r>l/m — - 



47r TT 



(17) 



Here the renormalized coupling constant is defined at //^ = 0. 

Comparing the theoretical expression (17) with the experimental data (5)-(6) allows to 
establish the numerical value of the expansion parameter a. In the first order of perturbation 
theory, i.e. retaining only e-term and neglecting e^-term in (17), we obtain 



g e ^ , Q2 



^o(r) = — > e ^ Q ^ a 



r r ' "'^-^ 47r"^""" in 137 

The value 

.2 1 



e 



47r"'"" 137 



< 1 (19) 



justifies assumption e < 1 used in deriving expression (18), when we have neglected term ~ e^. 
If more accuracy is desired, then e^-term is retained and the value of e emerges as a solution 
of a cubic equation. 

After the value of the expansion parameter is established, it is possible to find its value for 
a different subtraction point /i. This is done by analyzing the renormalization group equation 
[5]: 

ae(^_e^ 

To solve equation (20), the boundary condition is necessary, and this boundary condition is not 
given by the renormalization group. This boundary condition, numeric value of renormalized 
coupling constant 

e\0) _ J_ 
47r ~ 137 

is established from the comparison of the theoretical expressions with the experimental data as 
it is described above. 



Discussion 

Let us summarize the procedure of establishing the numeric value QED coupling constant. 
Prescription can be described as follows: 

• Choose physical quantity a. Numeric value of measurable a is known from the experiment. 

• Calculate a in the framework of renormalized Quantum Field Theory. In QED with the 
two input parameters the expression is (7(m, e)th- 

• Provided that the value of m is known, find e from the equation 

(T(m, e)th = o- 
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For (7 we have chosen Q, the coefficient of 1/r appearing in the Coulomb law. As it is mentioned 
earlier, any physical quantity will do fine. E.g., one can determine numeric value of e from the 
measurement of anomalous magnetic moment of electron [7] or quantized Hall conductivity 
[8]. Choosing Coulomb law seems bit more intuitive and transparent since there is no need 
to involve measurable depending on purely quantum features of electrodynamics (anomalous 
magnetic moment) or to apply QED to a problems of mesoscopic physics (Quantum Hall Effect). 
Of course, any choice for a yields the same result as in (19), just the accuracy may be different. 
Establishing the value of the QED coupling constant with the best possible accuracy was not 
our aim here - this goal requires some "heavy-weight" calculations, e.g., calculation of the 
anomalous magnetic moment of the muon in higher orders of perturbation theory, accounting 
for all the degrees of freedom of the Standard Model (electron, muon and tau leptons, quark, 
and W and Z bosons) [9]. 
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